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Abstract 

We construct a simple two-phase model of the nucleon structure 
functions valid for both small and large Q 2 and in the broad range 
of Bjorken x. The model incorporates hadron dominance at small 
x and Q 2 and parton model at large Q 2 . The VDM contribution is 
modified for small fluctuation times of the hadronic state of the photon. 
With two free parameters we describe SLAC, CERN NMC, Fermilab 
E665 and CERN BCDMS data for both proton and deuteron structure 
functions. Our model explains some phenomenological higher- twist 
effects extracted from earlier analyses. A good description of the NMC 
F^{x) — F£{x) data is obtained in contrast to other models in the 
literature. We predict faster vanishing of the partonic component at 
low Q 2 than previously expected and strong Q 2 dependence of the 
Gottfried Sum Rule below Q 2 « 4 GeV 2 . 

1 Introduction 



The standard deep inelastic scattering picture applies when the four- 
momentum transfer squared from the lepton line to the hadron line (Q 2 ) 
is large. When virtual photon wave length increases and reaches the size 
of the nucleon one may expect a transition to another regime where the 
standard partonic model is no longer valid. In this region a kinematical 
constraint Q guarantees the vanishing of the F2(x,Q 2 ) structure function. 
This requirement is not embodied in the perturbative parton distributions. 
A phenomenological fit based on parton screening was proposed in M to 
satisfy this condition by introducing an extra form factor. 
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The recent low-Q 2 data from HERA Q have triggered many phe- 
nomenological analyses. Especially interesting is the unexplored transition 
region. At present there is no consensus on the details of the transition 
mechanism. Here we concentrate on the region of somewhat larger x rather 
than the new HERA data. We shall demonstrate that also at larger x a 
similar transition due to vanishing partonic components at small Q 2 takes 
place, although it is not directly seen from experimental data. 

It is common wisdom that the vector dominance model applies at low 
Q 2 while the parton model describes the region of large Q 2 , leading at 
lowest order to Bjorken scaling, and to logarithmic scaling violation in higher 
orders of QCD. A proposal was made in Ref.(?J to unify both the limits in 
a consistent dispersion method approach. In the traditional formulation of 
the VDM one is limited to large lifetimes of hadronic fluctuations of the 
virtual photon, i.e. small Bjorken x < 0.1 for the existing data. It is a 
purpose of this paper to generalize the model to a full range of Q 2 and x 
by introducing extra phenomenological form factors to be adjusted to the 
experimental data. 

Some authors believe that it is old-fashioned to talk about VDM con- 
tribution in the QCD era. However, VDM effects appear naturally in the 
time-like region in the production of vector mesons in e + e~ collisions. These 
effects cannot be described in terms of perturbative QCD, as in the produc- 
tion of resonances many complicated nonperturbative effects take place. The 
physics must be similar in the space-like region. We shall demonstrate that 
it is essential to include this contribution explicitly in order to describe the 
structure functions at low Q 2 . 

In the next section we outline our model and discuss how to choose its 
basic parameters in a model independent way. In section 3 we discuss a fit 
of the remaining parameters to the fixed target data and present results of 
the fitting procedure for the proton and deuteron structure functions. In 
addition we compare the result of our model for i 7 ^ — and some subtle 
isovector higher-twist effects with another low-Q 2 model. Finally we discuss 
some interesting predictions which could be verified in the future. 

2 The model 

As in Ref. jjj the total nucleon structure function is represented as a 
sum of the standard vector dominance part, important at small Q 2 and/or 
small Bjorken x, and partonic (part) piece which dominates over the vector 
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dominance (VDM) part at large Q 2 : 

F 2 N (x, Q 2 ) = F?> VDM (x, Q 2 ) + F 2 N ' part (x, Q 2 ) . (1) 

The standard range of applicability of vector dominance contribution is lim- 
ited to large invariant masses of the hadronic system (W), i.e. small values 
of x. In the target (nucleon) reference frame the time of life of the hadronic 
fluctuation is given according to the uncertainty principle as r ~ 1 /AE with 



AE = .JM 2 + |q|2 - + | q |2 ( (2 ) 

where My is the mass of the hadronic fluctuation (vector meson mass). In 
terms of the photon virtuality and Bjorken x this can be expressed as 

AE = Ml±9l . Mn x . (3) 

Q 

As the energy transfer v — > oo the time of life of the hadronic fluctuaction 
becomes r ~ -r-r^- . It is natural to expect small VDM contribution when 

the time of life of the hadronic fluctuation is small. We shall model this fact 
by introducing a form factor £1(t) = Q(x,Q 2 ). Then the modified vector 
dominance contribution can be written as: 

2 [ >~ vr ^ 1 2 V {Q 2 + M 2 ) 2 Uv ^^>- W 

In the present paper we take 7's calculated from the leptonic decays of vector 
mesons which include finite width corrections || 7?/47r = 2.54, 'y 2 /^ = 
20.5, 7 2 /4tt = 11.7. Q 

In general one can try different functional forms for S7. In the present 
analysis we shall use only exponential and Gaussian form factors 

Q(x,Q 2 ) =exp(-AE/X E ) , 
n(x,Q 2 ) = exp(-(A£/A G ) 2 ) . (5) 



As in Ref.0 we take the partonic contribution as 

Q 2 

Q 2 + Ql 



F?> part (x, Q 2 ) = ®\ 2 ■ F« symp (x, Q 2 ) . (6) 



where x = w2 _Z^2 +Q , and Q 2 = Q 2 + Q 2 . The F^ symp (x,Q 2 ) above 



"N 



denotes the standard partonic structure function which in the leading order 
can be expressed in terms of the quark distributions: 

1 Please note different normalization of 7's in comparison to In. 
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Fr mp (x, Q 2 )=x- Ef ej ■ [q f (x, Q 2 ) + q f (x, Q 2 )] . 

The extra factor in front of Eq.@ assures a correct kinematic beheviour in 
the limit Q 2 — > 0. In general Qq, Q 2 and Q\ can be slightly different. In 
the following section we shall consider different options. 

At large Bjorken x one has to include also the so-called target mass cor- 



rections. Their origin is mainly kinematic [1C]. In our approximate treat- 
ment we substitute the Bjorken variable x in the partonic distributions by 
the Nachtmann variable £ jll| given by: 



2r 

1 + a/(1 + 



which is the dominant modification. 

In principle i^ S2/mp (x, Q 2 ) could be obtained in any realistic model of the 
nucleon combined with QCD evolution. We leave the rather difficult problem 
of modeling the partonic distributions for future studies. We expect that at 
not too small x > 0.01, the region of the interest of the present paper, the 
leading order Gliick-Reya-Vogt (GRV) parametrization of F2' asymp (x,Q 2 ) 
and F2' asymp (x,Q 2 ) should be adequate. Furthermore in our opinion the 
parametrization j|] with the valence-like input for the sea quark distribu- 
tions and d - u asymmetry built in incorporates in a phenomenological way 
nonperturbative effects caused by the meson cloud in the nucleon jl2| . 

The total cross section for (vector meson) - (nucleon) collision is not 
well known. Above meson-nucleon resonances, one may expect the following 
approximation to hold: 



—tot I —tot 

a 7T+ P + a TT-p 



tot tot _ ' 

a pN — a ujN — 2 

—tot ~tot I —tot ' 

a <t>N - a K+p + a K-p 



„tot i —tot 

a TT+ P "+* a %-p 



(8) 



Using a simple Regge-inspired parametrizations by Donnachie-Landshoff 
[13 1 of the total ttN and KN cross sections we get simple and economic 



parametrizations for energies above nucleon resonances 

af N = tffif = 13.63 • s - 0808 + 31.79 • s" ' 4525 , 

af N = 10.01 . s °-°8°8 + 2.72 . s -°- 4525 , (9) 

where the resulting cross sections are in mb. 
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We expect that our model should be valid in a broad range of x and Q 2 
except for very small x < 0.001, where genuine effects of BFKL pomeron 
physics could show up, and except for very large x, where the energy (s 1 / 2 in 
Eq.@) is small and the behaviour of the total VN cross section is essentially 
unknown. Because our main interest is in the transition region, the large Q 2 
data were not taken into account in the fit. There the partonic contribution 
is by far dominant and the GRV parametrization H is known to provide a 
reliable description of the data. 



3 Comparison to experimental data 

Most of the previous parametrizations in the literature centered on the 
proton structure function. In the present analysis we are equally interested 
in both proton and neutron structure functions. Achieving this goal requires 
special selection of the experimental data with similar statistics and similar 
range in (x,Q 2 ) for proton and deuteron structure function. In Fig.l we dis- 
play the experimental data for proton (left panel) and deuteron (right panel) 
structure functions chosen in our fit. We have selected only NMC, E665 and 



SLAC sets of data [18] for both proton and deuteron structure functions, 
together amounting to 1833 experimental points: 901 for the proton struc- 
ture function and 932 for the deuteron structure function. According to 
the arguments presented above, we have omitted BCDMS and HERA data 
in the fitting procedure but these will be compared to our parametrization 
when discussing the quality of the fit. 

The deuteron structure function has been calculated as 

Fi (x, Q 2 ) = \ M[x, Q 2 ) + F?(x, Q 2 )] , (10) 

i.e. we have neglected all nuclear effects such as shadowing, antishadow- 
ing due to excess mesons, Fermi motion, binding, etc, which are known to 
be relatively small for the structure function of the deuteron [|l3, [15], [ill], 
one of the most loosely bound nuclear systems. In addition we have as- 
sumed isospin symmetry for the proton and neutron quark distributions, 
i.e. u n (x, Q 2 ) = d p (x, Q 2 ), d n (x, Q 2 ) = u p (x, Q 2 ) and s n (x, Q 2 ) = s p (x, Q 2 ). 
The charm contribution, which in the GRV parametrization || is due to the 
photon-gluon fusion, is in practice negligible in the region of x and Q 2 taken 
in the fit, and therefore is omitted throughout the present analysis. 

The results of the fit are summarised in Table 1. Because in general 
Qq, Q\ and Q\ can be different, there are 4 independent free parameters 
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of the model. In order to limit the number of parameters we have imposed 
extra conditions as specified in Table 1. A series of seven fits has been 
performed. In all cases considered the number of free parameters has been 
reduced to two: the cut-off parameter of the form factor and Qq. Statistical 
and systematic errors were added in quadrature when calculating x 2 ■ Only 
data with Q 2 > 0.25 GeV 2 were taken in the fit which is connected with 
the domain of applicability of the GRV parametrization. The values of the 
parameters found are given in each case in parentheses below the value of x 2 
per degree of freedom. In addition to the combined fit, which includes both 
proton F% and deuteron i 7 ^ structure function data, we show the result of 
the fit separately for proton and deuteron structure functions. As can be 
seen from the table fairly similar values of parameters are found for the 
proton and deuteron structure function and the x 2 P er degree of freedom is 
slightly worse in the latter case which can be due to the omission of nuclear 
effects as mentioned above. The best fit (No 1 in the table) is obtained 
with Qi = Q 2 . = (fits of similar quality can be obtained with very small 
values of Q\ ~ 0.1 GeV 2 and Q\ ~ 0.1 GeV 2 ). While the value of X 2 
does not practically depend on the type of the form factor (exponential vs. 
Gaussian), much larger value of Qq is found for the Gaussian (Qq = 0.84 
GeV 2 ) than for the exponential (Qq = 0.52 GeV 2 ) parametrization. The 
value of Qq found here is smaller than in the original Badelek-Kwiecihski 
model 0] but larger than that found by HI collaboration in the fit to low-x 
data [[|. 

Although the resulting x 2 is similar in both cases, the F^ix^/F^ (x) 
ratio for x — > 1 prefers the Gaussian form factor. While the vector meson 
contribution with the exponential form factor survives up to relatively large 
x, with the Gaussian form factor it is negligible at large x.. 

For comparison the GRV parametrization of quark distributions alone 
yields: 

X 2 ' /Ndof = 9.74 (21.48) (proton structure functions), 
X 2 ' /Ndof = 13.73 (32.99) (deuteron structure functions), 
x 2 /N dof = 11.77 (27.33) (combined data), 

where the first numbers include target mass corrections and for illustration 
in parentheses their counterparts without target mass corrections are given. 
Clearly the inclusion of the target mass effects is essential and only such 
results will be discussed in the course of the present paper. 

The agreement of the CKMT parametrization is comparable to that 
obtained in our model. For instance for parametrization (b) in Table 2 in 
the second paper ||, which includes new HERA data: 
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X 2 /Ndof = 2.22 (1.00) (proton structure functions), 
X 2 /Ndof = 3.54 (3.59) (deuteron structure functions), 
x 2 /N dof = 2.89 (2.33) (combined data), 

where in the parentheses we present x 2 f° r Q 2 < 4 GeV 2 i.e. in the region 
of applicability of the CKMT parametrization. ^] We note much better 
description of the proton data in comparison to the deuteron data. The 
agreement of the Donnachie-Landshoff parametrization || with the proton 
structure function data is of similar quality. 

In Fig. 2 we present for completeness a map of x 2 f° r our best fit as a 
function of model parameters Qq and A. A well defined minimum of % 2 
for Xq ~ 0.5 GeV and Qq ~ 0.85 GeV 2 can be seen. The experimental 
statistical uncertainty of the obtained parameters Xg and Qq is less than 1 
%. 

Some examples of the fit quality can be seen in Fig. 3 (x-dependence for 
different values of Q 2 = 0.585, 1.1, 2.0, 3.5 GeV 2 ) and in Figs.4, 5 {Q 2 - 
dependence for different values of Bjorken x = 0.00127, 0.0125, 0.05, 0.10, 
0.18, 0.35, 0.55, 0.75). Shown are experimental data 18] which differ from 
the nominal Q 2 or Bjorken x specified in Fig. 3, 4, 5 by less than ± 3 %. An 
excellent fit is obtained for Q 2 > 4 GeV 2 (not shown in Fig. 3), although the 
VDM contribution stays large up to 10 GeV 2 . In comparison to the GRV 
parametrization (dashed line) our model describes much better the region 
of small Q 2 < 3 GeV 2 , especially at intermediate Bjorken x: 0.05 < x < 0.4. 
The CKMT model (long-dashed line), shown according to the philosophy in 
H for Q 2 < 10 GeV 2 gives a better fit at very small Bjorken x. However, 
one may expect here a few more effects which will be discussed below. ^ It 
is however slightly worse as far as isovector quantities are considered, as will 
be discussed later. There seems to be a systematically small (up to about 5 
%) discrepancy between our model and the data for Q 2 < 2 GeV 2 and x = 
0.1 - 0.3. This is caused by some higher-twist effects due to the production 
of the 7riV |l7j and 7rA exclusive channels and will be discussed elsewhere. 
A fit of similar quality is obtained in our model for the proton (left panels) 
and deuteron (right panels) structure functions. Rather good agreement of 
our model with the BCDMS data can be observed in Figs.4 and 5 in spite 
of the fact that the data were not used in the fitting procedure. 

At very small x < 0.01 the description of the data becomes worse. This 

2 The number of experimental points is reduced then to 354 and 373 for proton and 
deuteron structure functions, respectively 

3 The use of the next-to- leading order structure functions || in our model would improve 
the description of \ow-x data, discussion of which is left for a separate analysis. 
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is partially due to the use of the leading order approximation. The fit to 
the fixed-target data prefers x ~ x and Q 2 sa Q 2 (see Table 1 and the 
discussion therein). On the other hand, the HERA data would prefer Q 2 7^ 
and Q\ 7^ 0. If we included the HERA data in the fit the description 
of the fixed target data would become worse. There are no fundamental 
reasons for the parameters in both regions to be identical. In addition at 
very small x other effects of isoscalar character, not included here, such as 



heavy long-lived fluctuations of the incoming photon [19| and/or BFKL 



pomeron effects [20], may become important. 

For illustration a VDM contribution modified by a form factor (|5|) is 
shown separately by the short-dashed line. The partonic component can be 
obtained as a difference between the solid and VDM line. It can be seen 
from Figs. 3-5 that the partonic component decreases towards small Q 2 . This 
decrease is faster than one could directly infer from the failure of the GRV 
parametrization at low Q 2 because in our model a part of the strength resides 
in the VDM contribution. The modified VDM contribution is sizeable for 
small values of Bjorken x and not too large Q 2 and survives up to relatively 
large Q 2 . At Q 2 > 3.5 GeV 2 the structure functions in our model almost 
coincide with those in the GRV parametrization despite that the VDM term 
is still not small. For Q 2 — > 00 only the partonic contribution survives and 
F 2 {x,Q 2 ) - Fr\x,Q 2 ) -» F?ny(x,Q 2 ). 

The deviations from the partonic model can be also studied in the lan- 
guage of higher-twist corrections. Then the structure function can be written 
formally as 



F^ n (x,Q 2 ) = F^ LT (x,Q 2 ) 



cg/"(x) <?j\x) , 



Q- 



Q 4 



Cii) 



However, in empirical analyses one usually includes only one term in (11) 



? p/n 



(x,Q 2 



-,p/n,LT 



(x,Q 2 



1 + 



n p/n 



(X) 



Q 2 



(12) 



In our model for Q 2 ~ My, Qq there are an infinite number of active terms 
in the expansion of the structure function ([II]) . Therefore the coefficient (?l n 
(the same is true for the deuteron counterpart c d ) in (^) becomes effectively 
Q 2 -dependent &l n (x) = c p / n {x,Q 2 ). 

As an example in Fig. 6 we show c p and c d as a function of Bjorken x 
for three different values of Q 2 = 2, 4, 8 GeV 2 in the range relevant for the 
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analysis in [23]. In order to correctly compare our results for c p and c d with 
the results of the analysis in |23| the structure function F^ n ' LT (x,Q 2 ) in 



Eq.(|12|) will include complete target mass corrections calculated according 
to Ref.[jlO|. A fairly similar pattern is obtained for c p and c d especially 
at small x. The rise of c p or c d for x — > 1 is caused by our treatment of 
the target mass corrections and partially by the VDM contribution which 
survives in our model up to relatively large x. We obtain small negative c p 
and c d for x < 0.3 in agreement with |23|]. The smallnest of c p and c d in our 
model for x < 0.3 is due to the cancellation of the positive VDM contribution 
and a negative contribution caused by the external damping factor q2 + q$ of 
the partonic contribution in Eq.(||). The CKMT parametrization (shown 
only in its applicability range for Q 2 = 2, 4 GeV 2 ) provides a very good 
explanation of c p . It predicts, however, somewhat larger negative c d for x < 
0.3. This will have some unwanted consequences for c p — c n discussed below. 

In Fig. 7 we show c p — c n for Q 2 = 2, 4 GeV 2 together with empirical 
results from |24j]. Rather strong Q 2 -dependence of c p — c n is observed. Our 
model correctly describes the trend of the experimental data. For compari- 
son we show also the result obtained by means of the CKMT parametrization 
(long-dashed lines) of the structure functions which somewhat fails to re- 
produce the details of the empirical results from j24|, especially for small 
Bjorken x. To our best knowledge no other model in the literature is able 
to describe quantitatively this subtle higher-twist effect. 

Our model seems to provide a very good description of some isovector 
quantities. As an example in Fig.8 we present F%(x, Q 2 ) — F%{x, Q 2 ) at Q 2 = 
4 GeV 2 obtained in our model (solid lines for different form factors), as well 
as the results obtained with the GRV parametrization (dashed line) and in 
the CKMT model (long-dashed line).Q The NMC data @] prefer rather our 
model. As a consequence of the imperfect description of the deuteron data 
the CKMT model fails to describe the difference F^{x) — F£{x) for x < 0.3. 
The success of our model is related to the violation of the Gottfried Sum 
Rule and /or d — u asymmetry which is included in our model explicitly. In 
contrast to our model in the CKMT model for Q 2 > 2 GeV 2 the Gottfried 
Sum Rule Sg = §■ 



No evolution of the CKMT quark distributions was done, but it is negligible between 
4 GeV 2 for the nonsinglet quantity. 
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4 Conclusions and discussion 



We have constructed a simple model incorporating nonperturbative struc- 
ture of the nucleon and photon. Our model is a generalization of the well 
known and successful Badelek-Kwiecihski model Q. While the original 
Badelek-Kwiecihski model is by construction limited to a small- x region, 
our model is intended to be valid in much broader range. The original VDM 
model assumes implicitly a large coherence length for the photon-hadron 
fluctuation, i.e. assumes that the hadronic fluctuation is formed far up- 
streem of the target. When the fluctuation length becomes small the VDM 
is expected to break-down. This effect has been modelled by introducing an 
extra form factor. As a result we have succeeded in constructing a physically 
motivated parametrization of both proton and deuteron structure functions. 
In comparison to the pure partonic models with QCD evolution our model 
leads to a much better agreement at low Q 2 in a broad range of x. 

With only two free parameters we have managed to describe well the 
transition from the high- to low-Q 2 region simultaneously for the proton 
and deuteron structure functions. Our analysis of the experimental data 
indicates that the QCD parton model begins to fail already at Q 2 as high 
as about 3 GeV 2 . This value is larger than commonly believed. 

In our discussion we have omitted the region of the HERA data. In our 
opinion the physics there may be slightly more complicated. Other effects 
of isoscalar character, not included in our analysis, for example the BFKL 
pomeron effects p0|, may become important. 

In contrast to other models in the literature we obtain a very good 
description of the NMC — F<£ data [21] where the previously mentioned 
isoscalar effects cancel. 

Recently an intriguing, although small, difference between d — u asym- 
metry obtained from recent E866 Drell-Yan data [22| and muon DIS NMC 
data [21] has been observed. At least part of the difference can be under- 
stood in our model. We expect for Q 2 smaller than about 4 GeV 2 an extra 
Q 2 dependence of some parton model sum rules. We predict a rather strong 
Q 2 effect for the integrand of the Gottfried Sum Rule where in the first 
approximation the VDM contribution cancels. 

Recently in the literature there has been sizeable activity towards a 
better understanding of higher-twist effects. Some were estimated within 
the operator product expansion, some in terms of the QCD sum rules. It is, 
however, rather difficult to predict the absolute normalization of the higher- 
twist effects. Our model leads to relatively large higher-twist contributions. 
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For some observables, like structure functions, they almost cancel. For other 
observables, like i*f — Kp, the cancellation is not so effective. Our model 
provides specific higher-twist effects not discussed to date in the literature. 
This will be a subject of a future separate analysis. 
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Table 1: A compilation of the results obtained from our fit. The x 2 P^ r 
degree of freedom are given in first lines whereas the pairs of numbers in 
second lines are the parameters (X (GeV), Ql (GeV 2 )) found in the fit. 





exponential 


Gaussian 


case 


Ff and F$ 


F p 

r 2 




Ff and F£ 




n 


i) Qi = Qi = o 


2.39 
(0.31,0.52) 


2.13 
(0.30,0.51) 


2.63 
(0.31,0.53) 


2.38 
(0.50,0.84) 


1.90 
(0.49,0.77) 


2.66 

(0.50,0.87) 


2a) Qi = Ql 


3.24 

(0.34,0.51) 


2.79 
(0.33,0.50) 


3.68 
(0.34,0.51) 


2.79 
(0.53,0.77) 


2.27 
(0.53,0.77) 


3.15 
(0.53,0.84) 


2b) Qi =0.5 GeV z 


3.20 
(0.34,0.53) 


2.77 
(0.34,0.52) 


3.61 
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Figure 1: The presentation in the (x,Q 2 ) plane of the experimental data 
included in the fit for proton (left panel) and deuteron (right panel) structure 
functions. For reference shown are lines with constant W = 2 GeV which 
conventionally separate resonance and DIS regions. 
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Figure 2: A map of the x 2 V er degree of freedom in the (X, Qq) space of 
the model parameters for combined fit to the proton and deuteron structure 
functions with the Gaussian form factor. 
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Figure 3: Comparison of the model results with experimental data for 
(l.h.s.) and (r.h.s.) as a function of Bjorken x for different values 
of Q 2 = 0.585, 1.1, 2.0, 3.5 GeV 2 . The solid line corresponds to our 
full model with the Gaussian form factor. We present also the modified 
VDM contribution (short- dashed) and for comparison also the result obtained 
with GRV parametrization (corrected for target mass effects) of quark 
distributions (dashed line) and that\(jf the CKMT model (long-dashed 
line) . 
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Figure 5: The same as in Fig. 4 for different values of x = 0.18, 0.35, 0.55, 
0.75. 
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Figure 6: The twist-four coefficients c p and c d defined by Eq. p2] ) obtained in 
our model for Q 2 = 2,4,8 GeV 2 . The coefficients obtained from the CKMT 
parametrization are shown by the long-dashed line for Q 2 = 2,4 GeV 2 . The 
experimental data are taken from ft23(] . 
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Figure 7: The difference of twist-four coefficients c p — c n . The data points 
are from Ref. The solid lines are the results of our model while the 

dashed lines are obtained with the CKMT-parametrization. In both cases the 
two curves are for Q 2 = 2,4 GeV 2 . 
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Figure 8: F%(x,Q 2 ) - F 2 n (x,Q 2 ) at Q 2 =4 GeV 2 compared to the NMC 
data. The upper solid line corresponds to our model with the exponential 
form factor, the lower solid line to our model with the Gaussian form factor, 
the dashed line to the GRV parametrization and the long-dashed line to the 
CKMT parametrization. 
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